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Critical temperature of superconductor/ferromagnet bilayers 
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Superconductor/ferromagnet bilayers are known to exhibit nontrivial dependence of the critical temper- 
ature Tc on the thickness d/ of the ferromagnetic layer. We develop a general method for investigation of 
Tc as a function of the bilayer's parameters. It is shown that interference of quasiparticles makes Tc{df) a 
nonmonotonic function. The results are in good agreement with experiment. Our method also applies to 
multilayered structures. 

PACS: 74.50.+r, 74.80.Dm, 75.30.Et 



Recently, much attention has been paid to properties 
of hybrid proximity systems containing superconductors 
(S) and ferromagnets (F); new physical phenomena were 
predicted and observed in these systems |[ ^, . One 
of the most striking effects in SF layered structures is 
highly nonmonotonic dependence of the critical temper- 
ature Tc of the system on the thickness d/ of the ferro- 
magnetic layers. Experiments exploring this nonmono- 
tonic behavior have been performed previously on SF 
multilayers such as Nb/Gd @, Nb/Fe [|, V/V-Fe §, 
and Pb/Fe [|| but the results (and, in particular, com- 
parison between the experiments and theories) were not 
conclusive. 

To perform reliable experimental measurements of 
Tc{df), it is essential to have df large compared to the 
interatomic distance; this situation can be achieved only 
in the limit of weak ferromagnets. Active experimen- 
tal investigations of SF bilayers and multilayers based 
on Cu-Ni dilute ferromagnetic alloys are carried out by 
the group of Ryazanov^ j9|. In SF bilayers, they ob- 
served highly nonmonotonic dependence Tc{df). While 
the reason for this effect in multilayers can be the O-tt 
transition in a bilayer system with a single super- 
conductor this mechanism is irrelevant and the cause 
of the effect is quasiparticle interference specific to SF 
structures. 

In the present paper, motivated by the ROP experi- 
ment Q we study theoretically the critical temperature 
of SF bilayers. Previous theoretical investigations of Tc 
in SF structures were concentrated on systems with thin 
or thick S(F) layers [compared to the coherence length of 
the superconductor (ferromagnet)]; with SF boundaries 
having very low or very high transparency; besides, the 

^) e-mail: fominov@landau.ac.ru 

^^Ryazanov, Oboznov, Prokof'ev et al. — hereafter referenced 
as ROP. 



exchange energy was often assumed to be much larger 
then the critical temperature ||, |, |, |ll|, |l|| . The 
parameters of the ROP experiment do not correspond 
to any of these limiting cases. In the present paper we 
develop an approach giving opportunity to investigate 
not only the limiting cases of parameters, but also the 
intermediate region. Using our method, we find differ- 
ent types of nonmonotonic behavior of Tc as a function 
oi df such as minimum of Tc and even reentrant super- 
conductivity I p^ . Comparison of our theoretical predic- 
tions with the experimental data shows good agreement. 

We assume that dirty-limit conditions are fulfilled, 
and calculate the critical temperature of the bilayer 
within the framework of the linearized Usadel equations 
for S and F layers (the domain Q < x < ds'is occupied by 
the S metal, —df < x <Q — by the F metal, see Fig.l). 
Near Tc the normal Green function is G = sgnw„, and 
the Usadel equations for the anomalous function F take 
the form: 



^s^Tc: 



d'F, 
dx'^ 



|w„|F, + A^0, 0<x<d,- (1) 



d^F 

ttTc,— ^ - {\uJn\+iE^^^gnuJn)Ff = 0, (2) 



dx 



-df < X < 0; 



A\n^ = ttTY ( - Fs 
T ^ \ w„ 



(3) 



where = ^DJ2ttTcs, = ^Df/2^Tcs, w„ = 
TiT{2n + 1) with n = 0, ±1, ±2, . . . are the Matsubara 
frequencies, E'ox is the exchange energy, and Tcs is the 
critical temperature of the S layer. denotes the 

function F in the S(F) region. 
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Equations must be supplemented with the 

boundary conditions at the outer surfaces of the bilayer: 



dFsjds) ^ dFf{-df) 
dx dx 

as well as at the SF boundary: 



dF,{Q) ^ 
dx 

dFfiO) 



dx ' 



076 



dx 



F,(0)-F/(0), 7b 



Ps£.s 

RbA 



(4) 

(5) 
(6) 



Here psj are the normal state resistivities of the S and 
F metals, Rb is the total resistance of the SF boundary, 
and A is its area. The Usadel equation in the F layer is 
readily solved: 



Ff = C(w„) cosh {kf{x + df)) , 
1_ j \uJn\_\_iEc^sgnuj,, 



with k 



f 



(7) 
(8) 



and the boundary condition at a; = can be written in 
closed form with respect to Fs'. 



dFsjO) 
dx 
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76 + Bf{uj„) 



(9) 



where i3/(w„) = [kf£^fta,iih{kfdf)]~^. 

This boundary condition is complex. In order to 
rewrite it in real form, we do the usual trick and go 
over to the functions F^ = F([x)„) ± F(— cj„). The sym- 
metric properties of F^ and F~ are trivial, so we will 
treat only positive w„. The self-consistency equation is 
expressed only via the symmetric function F^: 



tj„>0 



(10) 



and the problem of determining Tc can be formulated 
in closed form with respect to F+. This is done as fol- 
lows. The Usadel equation for Ff does not contain A, 
hence it can be solved analytically. After that we ex- 
clude Ff from the boundary condition ^ and arrive 
at the effective boundary conditions for F+: 



^^^Ein = w{c.„)F+iO), 



dx 



^5it) = o, (11) 

ax 



where 



= 7- 



A,{-fb + 'ReBf) + -f 
As\jb + Bf\'^ + j{jb + ReBf) 

As ^ ksS^stanh(ksds), with fc. 



(12) 



1 luj^ 



£,s V ""^c; 



The self-consistency equation (p^), the boundary con- 
ditions ([Tll)-([l^) together with the Usadel equation for 
F+: 



d'F+ 
dx"^ 



2A = 



(13) 



will be used below to find the critical temperature of 
the bilayer. 

The Green function (in mathematical sense) of the 
problem ([ri] )-(p^ can be expressed via solutions wi, V2 
of Eq. (|3|) without A, satisfying the boundary condi- 
tions at a; = and x — ds, respectively: 



G{x,y;ujn) = - 



ks £,s/ 



saih.(ksds) + a ■ cosh (kgds) 

I vi{x)v2{y),x < y 



i'2{x)vi{y),y < x 



(14) 



where a — iy(w„)/fcs^s and 



Vi{x) ~ cosh{kgx) + a ■ smh{kgx), 
V2ix) = cosh {ks{x — ds)) . 



(15) 



Having found G{x,y;LOn), we can write the solution of 
Eqs. (©-(Ill) as 



F+{x;uJn) 



G(x,y;^„)A(y)dy. (16) 



Substituting this into the self-consistency equation 
we obtain 



A(x) In ^ 



= 27tT, J2 

uj„>0 L 



A{x) 



G{x,y;uJn)A{y)dy 



(17) 



This equation can be expressed in symbolic form: 
A\n{Tcs/Tc) = LA. Then is determined from the 
condition 



det L - 1 In ■ 







(18) 



that Eq. ( p7| ) has a nontrivial solution with respect to 
A. Numerically, we put our problem ([l7|)-(p^ on the 
spatial grid so that the linear operator L becomes a fi- 
nite matrix. 

Equations (p^)-(p^) are our central result; substitut- 
ing the concrete parameters of the system we can easily 
find the critical temperature numerically and in certain 
cases analytically. (The models considered previously 
|, [?[ § |ll, |ll|, |l|, |l|l correspond to the limiting cases 
of our theory.) 
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We apply our method to fit the ROP experimental 
data 1^; the result is presented in Fig. 2. Estimating 
the parameters ds = 11 nm, Tcs = 7K, = 7.5fiflcm, 

= 8.9 nm, Pf = 60 /ill cm, = 7.6 nm, 7 — 0.15 
from the experiment and fitting only £^ex and 76, we 
find good agreement between our theoretical predictions 
and the experimental data. The fitting procedure was 
the following: first, we determine i?cx ~ 130 K from 
the position of the minimum of Tc(d/); second, we find 
7b « 0.3 from fitting the vertical position of the curve. 
The deviation of our curve from the experimental points 
is small; it is most pronounced in the region of small df 
corresponding to the initial decrease of Tc. This is not 
unexpected because when df is of the order of a few 
nanometers, the thickness of the F film may vary signif- 
icantly along the film (which is not taken into account 
in our theory) and the thinnest films can even be formed 
by an array of islands rather then by continuous mate- 
rial. At the same time, we note that the minimum of T^. 
takes place at « 5nm, when with good accuracy the 
F layer has uniform thickness. 

The position of the minimum of Tc{df) can be esti- 
mated from qualitative arguments based on interference 
of quasiparticles in the ferromagnet. Let us consider a 
point X inside the F layer. According to Feynman's in- 
terpretation of quantum mechanics [l5[| , the quasiparti- 
cle wave function [we are interested in anomalous wave 
function of correlated quasiparticles, which character- 
izes the superconductivity; this function is equivalent 
to the anomalous Green function F{x)] may be repre- 
sented as a sum of the wave amplitudes over all classical 
trajectories; the wave amplitude for a given trajectory 
equals exp(iS'), where 5* is the classical action along this 
trajectory. To obtain our anomalous wave function we 
must sum over trajectories that (i) start and end at the 
point X, (ii) change the type of the quasiparticle (i.e., 
convert an electron into a hole or vice versa) . There are 
four kinds of trajectories which should be taken into ac- 
count (see Fig.l). Two of them (denoted 1 and 2) start 
in the direction toward the SF interface (as an electron 
and as a hole), experience the Andreev reflection and 
return to the point x. The other two trajectories (de- 
noted 3 and 4) start in the direction away from the in- 
terface, experience normal reflection at the outer surface 
of the F layer, move toward the SF interface, experience 
the Andreev reflection there, and finally return to the 
point X. The main contribution is given by the trajecto- 
ries normal to the interface. The corresponding actions 
are = —S2 = —Qx and S3 — —S4 ~ —Q{2df + x) 
(note that 2: < 0), where Q is the difference between the 
wave numbers of the electron and the hole. To make 
our arguments more clear, we assume that the ferro- 



magnet is strong, the SF interface is ideal, and con- 
sider the clean limit first: in this case Q — — kh = 
y/2m{E + E^^ + fi) - y/2m{-E - E^^ + /x) « 2E^^/v, 
where E is the quasiparticle energy, /i is the Fermi en- 
ergy, and V is the Fermi velocity. Thus the anoma- 
lous wave function of the quasiparticles is F{x) oc 
exp(iS'„) (X cos{Qdf) cos{Q{df + x)). The sup- 
pression of Tc by the ferromagnet is determined by the 
value of the wave function at the SF interface: F{Q) (x 
cos^{Qdf). The minimum of Tc corresponds to the min- 
imal value of F{Q) which is achieved at df = Tr/2Q. In 
the dirty limit the above expression for Q is replaced 
by Q = ■\/ Ecx/ Df, hence the minimum of Tc{df) takes 
place at 

In the case of the ROP bilayer |^ we obtain c?^'"'"^ « 
7nm, whereas the experimental value is 5nm (Fig. 2); 
thus our qualitative estimate appears to be reasonable. 

The method developed in this paper applies directly 
to multilayered SF structures (in particular, to trilayers) 
in the 0-state, where an SF bilayer can be considered as 
an elementary cell of the system. A generalization can 
be made, which allows to take account of possible su- 
perconductive and/or magnetic 7r-states. 

In conclusion, we have developed a method for cal- 
culating the critical temperature of a SF bilayer as a 
function of parameters of the junction. The approach 
developed here gives an opportunity to evaluate Tc in 
wide range of parameters. We demonstrate that there 
is good agreement between the experimental data and 
our theoretical predictions. Qualitative arguments are 
given, which explain the nonmonotonic behavior of the 
function Tc{df). Extensive details of our study will be 
published elsewhere Q. 
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Fig.l. The SF bilayer. The F and S layers occupy 
the regions —df<x<0 and < x < dg, respec- 
tively. The four types of trajectories contributing (in 
Feynman path integral sense) to the anomalous wave 
function of correlated quasiparticles are shown in the 
ferromagnetic region. The solid lines correspond to elec- 
trons, the dashed lines — to holes; the arrows indicate 
the direction of the velocity. 

Fig. 2. Theoretical fit to ROP's experimental data 
. In the experiment, Nb was the superconductor (with 
ds = 11 nm, Tcs = 7K) and Cuo.43Nio.57 was the weak 
ferromagnet. From our fit we estimate E'ox « 130K and 
76 « 0.3. 



